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Abstract 

We consider the case of bulk photons in a Lorentz violating brane background, 
with an asymmetric warping between space and time warp factors. A perturbative 
analysis, in a previous work, gave an energy dependent phase (or group) velocity of 
light: Vph{io) = Vph{0) — Cq ^'^ (Cg > 0), which was derived up to second order of 
time independent perturbation theory. In this paper, we go beyond the perturbative 
result and we study the nonperturbative behavior of the phase velocity for larger 
energies, by solving numerically an eigenvalue problem for the wave function of the 
zero mode (4D photon). In particular we see that V^/i(w) is in general a monotonically 
decreasing function which tends asymptotically to a final value Vph{oo). We compare 
with the results of perturbation theory and we obtain a very good agreement in the 
range of small energies. We also present a wave function analysis and we see that in 
the nonperturbative sector of the theory (very high energies), the zero mode and the 
massive KK modes tend to decouple from matter localized on the TeV brane. 
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1 Introduction 

Lorentz symmetry is assumed to be an exact symmetry of nature. However, there are many 
exotic theories, mainly quantum gravity and string models, which predict Lorentz violating 
effects in the high energy limit, see [T| and references therein. Such an effect is an energy 
dependent velocity of photons, as in the context of quantum gravity models of space time 
foam, where the vacuum behaves like a medium with a nontrivial subluminous refractive 
index. Accordingly, different times of arrival are expected when photons with very high 
energies, which are emitted simultaneously from remote astrophysical sources, reach the de- 
tectors of current experiments. Note that there are recent experimental data of MAGIC 
and FERMI telescopes, which imply a time delay of more energetic photons in comparison 
against lower-energy ones. However, such a difference may also have a conventional astro- 
physical interpretation: for example, photons with different energies may be emitted not 
simultaneously at their sources. 

Beyond quantum gravity models, an alternative mechanism which can produce a nontriv- 
ial vacuum refractive index based on brane models with an asymmetric space-time warping, 
was proposed. This mechanism was studied in [2j using time independent perturbation 
theory. Here we will study this mechanism in the nonperturbative regime, by solving numer- 
ically the eigenvalue problem for the wave function of the 4D photon. Such a study is useful 
for an analysis of extremely high energy cosmic phenomena in the context of asymmetric 
warp models. An example might be ultra high energy photons with energies higher than 
10^9 eV. 

Brane world models [31 HJ E] are models with extra dimensions which are used by theorists 
in order to address the hierarchy problem. According to this scenario, standard model 
particles are assumed to be localized in a three dimensional brane (our world), while gravitons 
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can propagate in the multidimensional bulk. Beyond the ADD scenario [41, |5j where the extra 
dimensions are assumed to be large, brane models in which the bulk space time is warped 
have been also proposed [6l[7]. In the case of warped space-time, the extra dimensions could 
be: (1) finite, if a second parallel brane world lies at a finite bulk distance from our world |6] 
or (2) infinite, if our world is viewed as an isolated brane, embedded in an (infinite) bulk 
space [7]. The previously mentioned model with the two branes, is often called first Randall 
Sundrum model (RSI-model). Generalizations of the above generic models, including, for 
instance, bulk fields along the extra dimension(s) or higher-order curvature corrections, have 
been also considered, see for example Refs. [H [9l [TOl |1T] and references therein. 
We will adopt the following generic ansatz for the metric in five dimensions 

ds^ = -a\z)de + f3^{z)dyi^ + -i^{z)dz^ (1) 

where z parameterizes the extra dimension. In contrast to the RS-model where the space 
and time warp factors are equal, in models with an asymmetric warping we have in general 
a{z) 7^ (3{z). Thus, although the induced metric on the brane (locahzed at 2 = for 
example) is Lorentz invariant upon considering the case a(0) = /3(0), the metric of Eq. ([T]) 
does not preserve 4D Lorentz invariance in the bulk since a{z) 7^ j3{z) for ^ 7^ 0. In such 
models Lorentz violation is due to bulk particles which can "feel" the difference between 
the space and time warp factors toward the extra dimension. In the standard brane-world 
scenario only gravitons are allowed to propagate in the bulk, hence, in the tree level, Lorentz 
violation effects are expected only in the gravitational sector. In Refs. [121 US] specific 
asymmetric models predict a superluminous propagation of gravitons. However, since the 
detection of gravitons is still not an experimental fact, we cannot use this effect in order to 
set restrictions to asymmetric brane models. 

As we have already mentioned, in Ref. [2] an asymmetric model where photons can 
freely move between two parallel branes in a 5D black hole background was considered. A 
perturbative analysis, of this model, gave an energy-dependent phase (or group) velocity of 
light: 

Vph{u:) = Vph{Q)-CGUJ^ (Cg>0) (2) 

and 

Vgr{uj) = Vgr{Q)-2CGU:'' (3) 

which was derived up to second order of time independent perturbation theory. Usually, in 
the conventional models the only bulk particles are the gravitons, but in the case of bulk 
photons which will be considered here we can set severe constraints to the free parameters 
of asymmetric models, see Ref. [2]. 

In this paper, we have examined the nonperturbative regime of the model in Ref. [2], 
by solving numerically an eigenvalue problem for the wave function of the zero mode (4D 
photon). We found that Vphij^) is indeed given by the perturbative formula of Eq. ([2]) in the 
range of small energies, but there is an inflexion point after which perturbation theory is not 
valid, then the phase velocity decreases monotonically with energy and tends asymptotically 
to a limiting value V^h(oo). 

In section 2 we introduce an asymmetric model with bulk photons, which consists of two 
branes in a 5D charged black hole background, and we examine in detail the corresponding 
junction conditions on the two branes. In section 3 we present the nonperturbative analysis 
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for the phase velocity and group velocity of 4D photon, while in section 4 we study the 
behavior of the wave function of the zero mode and the first KK excitation. Finally in 
section 5, we present our conclusions and we discus our results in connection with the MAGIC 
experiment, and the recent severe restrictions of ultra high energy cosmic rays on quadratic 
dispersion relations for the velocity of light. 



2 Asymmetric two brane models 

2.1 5D AdS-Reissner-Nordstrom black holes 

We consider an action which includes 5D gravity, a negative cosmological constant A, plus 
a bulk U(l) gauge field 



S 



MN 



B 



MN 



+ / d'x^^Cl:^^ + +J d'x^^p-\ 



(4) 

where is the five dimensional Newton constant, and Bmn = duHN — d^Hpi is the field 
strength of the U(l) gauge field Hm, with M, = 0, 1 ... 4. The four-dimensional terms in 
the action correspond to matter fields localized on the two branes of the model, which are 
located at r = r_|_ and r = r_ (r_ < r+), and described by two perfect fluids, localized on 
the two branes, with energy momentum tensors 
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As we will see in section 2.2, these brane terms are necessary for the solution of Eqs. 
and flTU]) below to satisfy the Israel junction conditions on the two branes. 
The corresponding Einstein equations can be written as 
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where the energy momentum tensor for the U(l) Gauge field is: 
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For the metric of the black hole solution we make the ansatz 

ds^ = -h{r)dt^ + r^r^dT? + h{r)-^dr^ 



(9) 



where dS^ = da"^ + a'^dQ'^ is the metric of the spatial 3-sections, which in our case are 



assumed to have zero curvature. Moreover, i is the AdS radius which is equal to y—j^- 
By solving the Einstein equations (j7]) we obtain: 



h{r) 



£2 



(10) 
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where /i is the mass (in units of the five dimensional Planck scale) and Q the charge of the 
5D AdS-Reissner- Nordstrom black hole. This, of course, presupposes the existence of extra 
bulk matter, namely a point-like source with mass /x and charge Q. Note that, in the case 
of nonzero charge Q, a non- vanishing component -Bor of the bulk field-strength tensor Bmn'- 

is necessary so that the solution satisfies the corresponding Einstein-Maxwell equations. 
2.2 Junction conditions 

We can use the 5D AdS-Reissner-Nordstrom black hole solution in order to construct two 
brane models. As a first step, we place two branes, one at the position r = r+ (Planck brane) 
and the other at the position r = r_ (TeV brane) (note that r+ > r_). We next assume that 
for r < r+ the 5D metric is given by Eq. iQ, while for r > r_|_ the metric is given by Eq. iQ 
upon the replacement r ^ r\/r. The metric which is obtained in this way is Z2-symmetric 
upon the replacement r ^ r\/r, and the points r_|_ and r_ correspond to the fixed points of 
the orbifold structure of the model. 

The next step is to glue the two independent slices of the metric by including two perfect 
fiuid energy momentum tensors on both branes, see Eqs. ([5]) and ([6]) above. Then we have to 
satisfy the junction conditions at the positions r = r+ and r = r_ (four junction conditions). 
In Refs. [T^ the junction condition for the corresponding single brane model has been 
derived. If we apply it in our case we take: 



6^/i(r±) = ±klp±r±, 18/i'(r±) = -kt{2 + 3u;±)p±r± (12) 
and after same algebra we obtain 

C - (l + i^1^Vi)^(l + ^W).* (13) 

If = {i + '-^ktfpi) = {i + '-^ktepi)e^ (14) 

where = SnG^, and (p+, p_) and (^4., p_) are the energy densities and pressures on the 
Planck and the TeV brane correspondingly . The equations of state are parameterized as 
usual by: 

^+=P+/p+, c^-=P-/p- (15) 
Note that p+ > (positive tension brane) and p_ < (negative tension brane). The 
parameter e is defined as: 

r_ 

6 = - (16) 
r+ 

In order to address the hierarchy problem, in a similar way with that of RSl-model, we have 
to choose a very large ratio e = 10~^^. Now, if we define the parameters 

P+ = P+, p- = ^p- (17) 
o o 

p£ _4 /=,2 Q ^ -6 /loN 

/^ = ^e , Q (18) 
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we obtain that 



/i = (l + c^+pt)e" =(l + ^-P-) (19) 

= ii+'-^PlK' = ii+'-^p'-) (20) 

In what follows we will consider that < p, « 1 and < << 1, as our purpose is to 
construct two brane models that are described by asymmetric metrics which are linearized 
perturbations around the RSI metric. See also Eq. (125]) for the perturbation Sh below. This 
implies that r+, which is the radius that determines the position of the Planck brane in the 
bulk, is (comparatively) a very large quantity 0. In particular, we have to satisfy both the 
inequalities r^e^ ^ ^/Jli and r^e^ ^ Qi. 

By solving Eqs. (fT9i) and (l20i) (they are four algebraic equations) we find the energy 
densities: 

pI = i-Sfie^ + 2Qh^ (21) 
pi = 1-3/2 + 2Q2 (22) 

and the equation of state parameters 

w+ = -1 - 2fx + 2Q^e^ (23) 
w_ = -1 - 2/i + 2Q2 (24) 

We see that the energy densities p+, p- and the state factor parameters w+, W- depend 
only on the constants p and and the hierarchy parameter e. Note that for = — 1 and 
= 1 we obtain the first RS-model. 

The equation of state parameters should respect the null energy condition > — 1, hence 
if we demand a;+ > — 1 and cj- > — 1 we obtain the following constraint : 

p < Q^e^ (25) 

This equation means that we can choose the parameters p (> 0) and (> 0) arbitrarily 
insofar as they satisfy the constraint of Eq. (!25|) . Note that between the two branes there 
are no horizons as the parameters p and are assumed to be very small, and the positions 
r+ and r_ of the two branes are very large. 

Finally, we would like to stress the fact that, although this problem has been considered 
previously in the literature, see Refs. [121 113], the results for the energy densities (Eqs. fl2Tl) 
and fl22l) ) and the state factors (Eqs. fl23|) and fl2^ ) are presented for the first time in the 
present paper. 



2.3 5D AdS-Reissner-Nordstrom Solution as a linearized pertur- 
bation around the Randall- Sundrum metric 

To write the 5D AdS-Reissner-Nordstrom solution as a linearized perturbation around the 
RS metric, we perform the following change of variables r — >■ z{r) in Eq. (Q: 

r = r+e~^ , for z > 

r = r+e^ ^ , for 2 < , (26) 
^Also r_ is a very large quantity because /I = |^ and — %Ji are assumed to be very small numbers. 
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If we rescale — >• ^a;^ (/i = 0, . . . , 3), we obtain: 

ds^ = -a\z)h{z)dt^ + a^{z)d^^ + h{z)-^dz'^ (27) 

where a{z) = e~*^'^', and k = is the inverse AdS^ radius. For the function h{z) we obtain: 

h{z) = l-5h{z), 5/1(2) = 3 e'^^l^l - 2 e^^l^l (28) 

The positions of the branes which are located at and r_ = r+e in the original coordinate 
system, are determined in the new coordinate system by the equations z = and z = Zc 
correspondingly, where e = e''^^" {tc = Zc/n is radius of the compact extra dimension). Note 
that the large hierarchy e ~ 10~^^ is achieved if we choose Zc — 37. In addition, we will 
assume that |5/i(2)| <^ 1 in the interval < z < Zc, or equivalently we adopt that 6h{z) is 
only a small perturbation around the RS-metric. We shall use the term Planck brane for 
the positive tension brane at the position z = 0, and the term TeV brane for the negative 
tension brane, at Zc- 



2.4 Bulk photons in asymmetric two brane models 

In this section we will study the case of a 5D massless U (1) gauge boson An in the background 
of an asymmetrically warped solution of the form of Eq. ( |271) . We stress that the gauge field 
An must not be confused with the gauge field H^, introduced in the previous section. As 
we will see later, we will identify the zero mode of An with the standard four dimensional 
photon. On the other hand Hn is an additional bulk field which does not interact with the 
charged particles on the brane. The equation of motion for An reads: 

-^Om {V99''''9'''Fns) = , (29) 

with Fns = QnAs — OsAn, and A^, 5* = 0, 1, ... 5. In the background metric of Eq. (|27|1 . 
Eq. ([29D gives: 

- d,ia\z)h{z)d,A^) - V^A, + j^d^.Aj =0, j = 1, 2, 3 , (30) 

where we have assumed the Coulomb gauge condition: 

V-1=0, Ao = 0, A, = 0, (31) 

which is suitable for the case of a Lorentz violating background. On setting in Eq. (!29|) : 

A,(a;,2) = e^P-^X,(^), P, = i-iu,p) (32) 

we obtain 

- {a\z)h{z)d.x} + {p'-^]x = (33) 

where for brevity we have dropped the index j from x- Note that the spectrum of Eq. (133!) 
is discrete, due to the orbifold boundary conditions [6], x'(0) = and x'(^c) = (where the 
prime denotes a z-derivative). 
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We would like review here briefly the spectrum in the case of RSI-model {6h = 0), which 
consist of a zero mode plus an inflnite tower of massive KK modes. It suffices to mention 
that the nonzero eigenvalues are: 

m(f) = x„ n = 1,2,3,... (34) 

where Xn are the roots of the zeroth order Bessel function Jo{xn) = 0. On adopting kzc ~ 12, 
which is the standard choice in order to connect electroweak (ew) and Planck scales (Mp = 
e^^^mew) in a RS framework [6j, one obtains that: 

m(f) ~ TeV , n = 1, 2, . . . . (35) 

The corresponding eigenfunctions are: 

Xo^ = T^, No = ^c (36) 

xf = Tr^^'Ji{^e^l, iV„ = ^Ji(x„) (37) 



Nn ' k " " V2k 
where the coefficients Nq, Nn are defined by the normalization condition: 

X^^\z)xt\z)dz = (38) 



3 Phase and group velocity of photons: a nonpertur- 
bative analysis and a comparison with perturbation 
theory 

If we introduce the dimensionless variable y = kz m Eq. fl55]) we obtain 
where 

h{y) = l-5e'y{^Ca-2e^y), 5 = Qh\ Ca = ^ (40) 

and the boundary conditions now read: 

X'(0) = 0, x'(^^c) = (41) 

Note that < < 1 if we demand the null energy condition to be satisfied, see section 
2.2 and Eq. fl2^ above. We have also checked that the value of has not a significant 
impact for our numerical analysis, hence in what follows we will assume that Ca = 1, which 
is the case where c<j+ = —1, see Eq. ( !23|) above. Now, it is convenient to introduce a new 
parameter 

A = g2 = 5 e-6 (42) 
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Figure 1: The phase velocity of the photon Vph as function of co/k for kzc = 2 and A = 
1, 0.5, 0.1, 0.01. We see that as A decreases the perturbative region of Vph spreads towards 
larger energies uj/k. 

which can help us to estimate where perturbation theory fails as an approximation for 
solving Eq. fl39|) . In particular, even for A << 1, we can trust perturbation theory only 
if the energy of the photon cu is relatively small. For larger energies cu the term □ to'^Sh 
{6h = 1 — /i^A^l)in Eq. fl33|) (or in Eq. fl39l) ) cannot be assumed small. In this case 
perturbation theory breaks down and a nonperturbative approach is necessary. 

The nonperturbative analysis of this section is reduced to an eigenvalue problem of the 
second order differential equation (l39l) with the boundary conditions (HTI) . We will consider 
that IpI is fixed and will try to determine the energy u in order to satisfy the boundary 
conditions of Eq. (14T!) . Note that there is an infinite tower of energies Un (n=0,l,2..) which 
are solutions of the above mentioned eigenvalue problem for given momentum |p|. The first 
eigenvalue Uq corresponds to the zero mode and the remaining eigenvalues Un {n ^ 0) to 
the massive KK excitations. Thus, we kept |p| fixed and we integrated numerically El, for 
a number of values oi uj {uj > |p|) separated by a relatively small constant step 5uj, Eq. 
fl39|) with the initial condition x(0) = 1, x'(0) = 0. We observed that as we increased a;, 

^We have used that h = 1 — 6h, hence the term ^ in Eq. ((551) can be written as + Sh) — uj"^ + uj'^Sh. 
The term ui'^Sh in the previous equation reveals the perturbative nature of Eq. [33] as i5/i « A <C 1. However, 
for large energies a; 1/ \/Sh perturbation theory breaks down. 

•^For the numerical solving of the eigenvalue problem we have used mathematica programming. 
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Figure 2: The phase velocity of the photon Vph as function of uj/k for kzc = 8 and A = 1,0.001. 
We see that even for larger hierarchies kzc (e = e~^^'') the qualitative features of the phase velocity 
do not change. 

with the constant step 6uj, the derivative x'{kzc) changed sign for first time (from positive 
to negative), in an interval of energies [uja,uja + Su]. In this interval the first eigenvalue ujq 
can be determined by bisection method. The infinite tower of energies Un {n ^ 0), can be 
determined in the same way. 
The phase velocity 

Vph = ^/IpI (43) 

of the photon (zero mode) as a function of u (measured in units of k) has been plotted in 
FigUl Fig 12] and Figj3l for several values of A and kzc (e = e'^^"). As we see, the phase 
velocity is a monotonically decreasing function which tends asymptotically to a constant 
value, seemingly equal to one independently from the parameters A and kz^. 
In Fig. IHwe have plotted the group and phase velocity as a function of u/k 



duj 
d\v\ 



Vgr = (44) 



We can observe that they have a very similar behavior. However, in contrast with the phase 
velocity, the group velocity becomes smaller than unity, then it increases and tends rapidly 
to unity. Also, as we see in Fig. HI the phase and the group velocity are equal in the low 
energy limit as it is expected. It is well known that the phase and group velocity in vacuum 
should be identical, however for larger energies we observe significant differences because 
the vacuum behaves as a medium with a non-trivial refractive index. The perturbative 
formulas for the phase and group velocity, of Eqs. ([2]) and ([3]), agree with our numerical 
analysis in the low energy limit, as they give the same value for the group and phase velocity 
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Figure 3: The phase velocity of the photon Vp^ as function olco/k for kzc = 3.5, 4, 5 and A = 0.01. 
As kzc increases the perturbative region shrinks in the low energy range u)/k « 1, and as a result 
the velocity Vph tends to its limiting value in a much faster way. 

for zero energy, decreasing quadratically with energy, while for larger energies outside the 
perturbative region they behave differently. 

Especially in Fig|5l we observe that there is an inflexion point Uf which separates the 
perturbative from the nonperturbative sector of the theory. In the perturbative sector {uo < 
ojf) we expect that 

Vph{i^) = Vph{Q) - Cg uj\ [Cg > 0) (45) 

This formula has been derived in Ref. ^ by using second order time independent pertur- 
bation theory. The parameters Vph{0) and Cq are given by the formulas of Eqs. (2.30) and 
(2.31) in [2], These formulas are suitable for numerical computations, if the parameters A 
and kzc are known. 

We also see, in FigEl that when uj crosses the inflexion point ujf the rate of decreasing of 
the phase velocity gets smaller and the phase velocity possesses an asymptotic value equal to 
one. Note also, that the perturbative range of the phase velocity increases (or the inflexion 
point u)f is displaced towards the right direction in the figures) for smaller values of the 
parameter A, as we see in FigH] and Figj2l 

At this point, we would like to stress that our numerical analysis is restricted to an 
unrealistic range of the parameter space of A and kZc, which is far beyond the physically 
interesting case of kzc — 37. However, it is reasonable that we are not in position to perform 
numerical computations in this it demands great accuracy. This is mainly due to 

the extremely large values of the exponential e^^^, that appear in Eq. fl25]) . 

Now, if we compare FigH] for kz^ = 2 and FigI2] for kzc = 8, we see that the qualitative 
features of the phase velocity as a function of u are unchanged, although there is a significant 
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Figure 4: The phase and the group velocity of the photon (Vph and Vgroup correspondingly) as a 
function of uj/k for kzc = 3 and A = 0.01. 

difference between the corresponding hierarchies ei = ~ 10^^ and £2 = ~ 10~^. We 
have also performed computations for larger kzc = 10 (£3 ~ 10~^) and we have confirmed 
the same behavior for the phase velocity. We also see that this behavior is independent from 
the parameter A, which determines the perturbative range of our model. Accordingly, in 
the conclusions we will consider an extrapolation assuming that the qualitative behavior of 
the phase velocity is also valid for e ~ 10~^^ which is the physically interesting case. 

Finally in Fig. [3] we see that as kzc increases the perturbative range of the phase velocity 
shrinks near the origin, where u/k <^ 1. This behavior is reasonable as uj in this figure is 
measured in units of k (or in units of Planck scale) and, in the case of realistic values of 
kzc — 37, the point where we have the breakdown of perturbation theory is expected to be 
several orders of magnitude smaller than the Planck scale (we will give an estimate of this 
point in conclusions). 

4 Wavefunction analysis 
4.1 Zero mode 

In Figl5] and Figj6] we have plotted the square of the normalized wave function measured 
in units of k {k'^Xnorm) several values of the energy u, assuming that the values of the 
parameters A and kzc are kept fixed. In particular, in FiglHl we observe that for small 
energies within the perturbative region of u {u < Uf) the wavefunction of the photon is 
almost constant. This is expected as in the case of RS-model {6h = or A = 0) the 
wavefunction can be obtained analytically and it is constant, see Eq. ( l36i) . As the energy of 
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Figure 5: The phase velocity of the photon Vph as function of uj/k for kzc = 3 and A = 0.01. The 
discrete points are the nonperturbative results, and the continuous line corresponds to perturbation 
theory: Vph{uj) = 1.00055 — 0.00012 {u>/k)'^, where the coefficients have been computed by the 
formulas of Eqs. (2.30) and (2.31) in Ref. [2]. We see that the inflexion point ujf ~ 1.1 k separates 
the perturbative from the nonperturbative sector of the theory. 



the photon increases the value of the wave function on the TeV brane {xnorm{zc)) decreases, 
see Figjni For lj = Uf the value of the wave function on the TeV brane is half of its value 
at ~ 0, while for larger values of (cu > ujf) it tends rapidly to zero. On the other 
hand, we see that the wave function on the Planck brane increases. Especially, in Fig. [7] 
we have plotted the wavefunction of the photon for even larger values of u, deep in the 
nonperturbative sector of the theory. It seems that the wavefunction on the Planck brane 
tends to take a limiting value, quite larger than that for ~ 0. 

If we take into account that the "effective" coupling constant of the zero mode (4D 
photon), with matter localized on the TeV brane, is proportional to Xnorm{kZc) (see p!5l[T6] ). 
we conclude that photons with very high energies u » ujf tend to decouple from matter 
which is localized on the TeV brane. Note, that a similar behavior has been observed for the 
massive KK modes of the 5D photon, which happens for even larger energies, as we see in 
the next section. On the other hand, in the case of Planck brane the coupling of the photon 
with matter increases with the energy, and tends to an asymptotic value. 

4.2 First KK excitation 

In this section we examine the wave function and the group velocity of the first KK mode. 
In particular, in Fig. [HI we see that the projection of the wave function on the TeV brane 
decreases as the energy of of the IKK mode increases, and for quite large uj the value of 
the normalized wave function becomes almost zero on the TeV brane. This means that 
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Figure 6: The square of the normahzed wave function k~^Xnorm of the zero mode (4D photon) as 
a function of kz for kzc = 3, A = 0.01 and to = 0.1k, 0.5k, 1.1k, 2k, 3k, 10k. We see that the value 
of the wave function on the TeV brane (which is proportional to the coupling of the zero mode with 
the localized matter on the TeV brane) tends to zero as the energy of the zero mode increases. On 
the other hand, on the Planck brane, we observe that the value of the wave function increases with 
the energy. 

for comparatively large energies where Lorentz violation effects become significant the IKK 
mode tends to decouple from matter, which is localized on TeV brane. Note, that a similar 
behavior was obtained for the zero mode in the previous section. In the case of higher KK 
excitations a similar behavior is expected. On the other hand the projection of the wave 
function on the Planck brane is almost constant independently from the energy u. 

The IKK mode is a massive particle, as for zero momentum p the energy u takes a 
nonnegative value u = ttiikk (t^ 0). In Fig. [HI for A = 0.01 and kzc = 3, we see that 
TT^iKK = 0.153/c, as it is the lower energy which is obtained for p = which corresponds 
to the inertial mass of the particle. Note that in the case of A = 0, where we can use the 
formula 

mi°) = x„fce-'=^% n = 1,2,3,... (46) 

that gives the masses of the KK excitations, where Xn are the roots of the zeroth order 
Bessel function Jo{xn) = 0. For n = 1 we obtain that mixK = 0.120/c. However we 
can use first order perturbation theory to correct this value (see ref. [2]), and now we 
obtain that rriiKK = 0.162/c which is close to the value rriiKK = 0.153/c that is obtained 
nonperturbatively. It is worth noting that in the realistic case kzc = 37 and A ~ 10~^ << 1 
these differences are expected to be more suppressed, and obviously are not detectable in 
the current high energy experiments, for example in LHC 

As the IKK mode is a massive particle the phase velocity is not suitable to describe its 
motion. For this reason in Fig, [9] we have plotted the group velocity of the particle as a 
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Figure 7: The square of the normahzed wave function k~^Xnorm of the zero mode (4D photon) as 
a function of kz for kzc = 3 and A = 0.01 and to = 5k, 15k, 25k, 35k. We observe that the value of 
the wave function, on the Planck brane, increases with the energy and tends to a constant value. 

function of u/k. We see that the energy has a lower bound which characterize the inertial 
mass of the particle as we have also explained in the previous paragraph. We also see that 
the group velocity is always smaller than unity, which is the standard velocity of light in the 
tree level of our model, and tends rapidly to this value Vgr = 1 as the energy u increases. 
Finally, we would like to note that in our model the group velocity of the zero mode, even if 
it becomes smaller than unity as we see in Fig. HJ is always larger than the group velocity 
of the first KK mode. 



5 Discussion 

We examined a two brane model where the 5D Lorentz invariance is spontaneously broken 
due to the nonstandard vacuum of a five- dimensional charged black hole. In this framework 
we found a mechanism which produces an energy dependent vacuum refractive index, as- 
suming that photons can freely move in the bulk, in contrast to the conventional brane world 
hypothesis. As perturbation theory was examined extensively in a previous work, in this 
paper we focused to the nonperturbative case by solving numerically the eigenvalue problem. 

We have mainly studied the phase (and the group) velocity of the zero mode, 4D photon, 
and we found that it is in general a monotonically decreasing function which for very large 
energies tends to unity, that is the standard velocity of light at tree level of our model. Note 
that a very similar behavior was obtained for the group velocity of photon, as we can see 
in Fig. m above. On the other hand, in the case of the first KK mode which is a massive 
particle, we found that the group velocity is always smaller than unity, and it cannot exceed 
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Figure 8: The square of the normahzed wave function k~^Xnorm IKK mode as a function of 

kz for kzc = 3 and A = 0.01 and uj = 0.5k, Ik, 2k, 3k, 5k, 10k. We see that the value of the wave 
function on the TeV brane tends to zero as the energy of the IKK-mode increases. On the other 
hand, on the Planck brane, we observe that the value of the wave function is almost constant. 

the group velocity of the zero mode in the high energy limit. 

By comparing with perturbation theory we found that there is an energy LOf after which 
perturbation theory breaks down. Specifically cu/ is the inflexion point (see Fig. [5]) where 
the quadratic dependence on energy terminates and the velocity tends to a limiting value. 
One could give an estimate of this point by comparing with the recent data of the current 
experiments of MAGIC H.E.S.S [12] and FERMI [20] telescopes. In the case of our 

model, which predicts a quadratic dependence on the energy for the velocity (see Eq. ([2])), 
the stringent bound is set by theMAGIC experiment: 

V = 1- i^—j , M2 > 2.6 X lO^^GeV (47) 

The above restriction was obtained in Ref. [H] by fitting the recent experimental data of 
MAGIC [T7] assuming a quadratic energy dependence for the photon refractive index. Hence, 
if we take the lowest bound for M2 (M2 = 2.6 x lO^^GeV) we conclude that the energy Uf, 
after which perturbation theory breaks down, should be quite smaller than 2.6 x lO^^Ge^. 
Note that the above upper limit, for the inflexion point Uf, if it is expressed in eV, gives a 
value equal to 2.6 x lO^^eV^ which is close (but smaller) to the energy range of the ultra high 
energy cosmic rays (particles with astrophysical origin and energies larger than the GZK 
limit 7 X 10^^ eV). Accordingly the quadratic dependence of velocity of light from the energy 
LJ, in the ultra high energy cosmic rays energy region, is not valid any more. Our analysis, 
in this region of energies, shows that the velocity of light is almost independent from the 
energy and it has taken a limiting value, which is the velocity of light at tree level of our 
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Figure 9: The group velocity Vgr{uj) as a function of energy in the case of the first KK mode, 
for kzc = 3 and A = 0.001. Note that the group velocity approaches the unity remaining always 
smaller than the group velocity of the zero mode. We also see in this figure that the energy range 
has a lower bound which is identified to the mass of the IKK mode. 

model. 

In Ref. [21] the authors found the following severe constraint for quadratic dispersion 
relations: 

V^=l-^6(^)', MpL = 10'' GeV, 6<2.4xlO-^ (48) 

which is due to the lack of observations of photons above the GZK limit, and appears to 
be several orders of magnitude stronger than the bounds of MAGIC observations (compare 
with Eq. fH7|) ). However, this constraint presupposes that the quadratic energy dependence 
of the velocity of light is valid in the ultra high energy cosmic ray energies, something which 
does not happen in our model as we have mentioned previously. We conclude that the 
constraints from the MAGIC observations are the most stringent ones which can be applied 
to our model. 
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